Existence of Temperature on the Nanoscale 
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We consider a regular chain of quantum particles with nearest neighbour interactions in a canonical 
state with temperature T. We analyse the conditions under which the state factors into a product of 
canonical density matrices with respect to groups of n particles each and under which these groups 
have the same temperature T. In quantum mechanics the minimum group size n m in depends on 
the temperature T, contrary to the classical case. We apply our analysis to a harmonic chain and 
find that n m in = const, for temperatures above the Debye temperature and n m i n oc T -3 below. 
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Recent progress in the synthesis and processing of 
materials with structures on nanometer length scales 
calls for better understanding of thermal properties of 
nanoscale devices, individual nanostracturcs and nanos- 
tructured materials. Experimental techniques have im- 
proved to such an extent that the measurement of ther- 
modynamic quantities like temperature with a spatial 
resolution on the nanometer scale seems within reach 
PI 0, 0- Q ■ These techniques have already been applied 
for a new type of scanning microscopy, using a tempera- 
ture sensor |fj that shows resolutions below lOOnm. 

To provide a basis for the interpretation of present day 
and future experiments in this held, it is indispensable 
to clarify the applicability of the concepts of thermody- 
namics to systems on small length scales. In this context, 
one question appears to be particularly important and in- 
teresting: Can temperature be meaningfully defined on 
nanometer length scales? 

The standard procedure to show that the thermody- 
namical limit and therefore temperature exists are based 
on the idea that, as the spatial extension increases, the 
surface of a region in space grows slower than its volume 
0. If the coupling potential is short-ranged enough, the 
interactions between one region and another become neg- 
ligible in the limit of infinite size. 

However, the full scaling behavior of these interactions 
with respect to the size of the parts has, to our knowl- 
edge, not been studied in detail yet 0, H, El 

For standard applications of thermodynamics this 
might not be very important since the number of par- 
ticles is typically so large that deviations from infinite 
systems may safely be neglected. Nevertheless, these dif- 
ferences could become significant as the considered sys- 
tems approach nanoscopic scales. Here it is of special 
interest to determine the "grain-size" needed to ensure 
the existence of local temperature. 

Since a quantum description becomes imperative at 
nanoscopic scales, the following approach appears to be 
reasonable: Consider a large quantum system, brought 
into a thermal state via interaction with its environment, 



divide this system into subgroups and analyse for what 
subgroup-size the concept of temperature is applicable. 

We adopt here the convention, that a local temper- 
ature exists, if the considered part of the system is in 
a canonical state, i.e. the distribution is an exponen- 
tially decaying function of energy characterised by one 
single parameter. This implies that there is a one to one 
mapping between temperature and the expectation val- 
ues of those observables, by which temperature is usually 
measured. After proper calibration, such measurements 
thus all yield the same temperature, contrary to distri- 
butions with several parameters. If the distribution were 
not exponentially decaying, physical quantities like en- 
ergy would not have "sharp" values. 

Our results arc thus obtained from calculations with 
canonical ensembles. Approaches based on different en- 
sembles (e.g. microcanonical ones), which are becing dis- 
cussed |l4f . may yield significantly different results. 

In this paper we study a chain of identical particles 
with next neighbour interactions which decomposes into 
Ng identical subgroups of n particles each. Assuming 
that the total system is in a thermal (canonical) state, 
we derive conditions on the interaction strength between 
the groups and the global temperature which ensure that 
the state of each group is approximately canonical. Fi- 
nally, we apply these conditions to a harmonic chain and 
calculate the minimal group size n m i ni which is found to 
depend on the global temperature and the desired accu- 
racy of the local thermodynamical description. 

We start by defining the Hamiltonian of our chain in 
the form, 



H 



with Hi 



■V(qi,q i+ i), (1) 



where the index i labels the particles of mass m and qi 
and pi are their positions and momenta, respectively. We 
assume periodic boundary conditions. 

We now form Ng groups of n subsystems each 
(index i —> (//— l)n + j; \x = 1, . . . , Ng; j = 1, ■ ■ ■ , n) and 



2 



split this Hamiltonian into two parts, 
H = Ha+I, 



(2) 



where Hq is the sum of the Hamiltonians of the isolated 
groups and / contains the interaction terms of each group 
with its neighbour group: 

JVo n 

Ho = / ] (Hfi — V(<lim, %n+\)) ; Hft = ^2^n(n-l)+j 

|U=1 3 = 1 
N G 

I = ^2V(q^ ni q^ n+1 ). (3) 

We label the eigenstates of the total Hamiltonian H and 
their energies with the greek indices (<p, tp) and eigen- 
states and energies of the group Hamiltonian H with 
latin indices (a, 6), 



H \<p) = E v \<p) and H \a) = E a \a). 



(4) 



Here, the states |o) are products of group eigenstates, 
with E a — Y^uSi Eft, where is the energy of one sub- 
group only, (Hfj, - Viq^n, q^n+i)) \a) = E M \a). 

We assume that the total system is in a thermal state 
with a density matrix, which reads 



(<p\p\ip) 



(5) 



in the eigenbasis of H . Here, Z is the partition sum and 
= (/cbT) -1 the inverse temperature. Transforming the 
density matrix © into the eigenbasis of Hq we obtain 



(o|p|a) = / w a (E) 

I Eg 



dE 



(6) 



for the diagonal elements. Here, the sum over all 
states | a) has been replaced by an integral over the en- 
ergy. The conditional probabilities w a (E) are given by 

w a (E v ) = r)(E v ) T,fHM=E v \ v )} l( a MI 2 with 7 l(E v ) de- 
noting the density of energy levels and the sum running 
over all states \<p) with H\<p) = E^cp). To compute the 
integral of equation (jSJ) we need to know the distribution 
of the conditional probabilities w a {E). 

In the limit of infinite number of groups Nq, w a (E v ) 
converges to a Gaussian normal distribution, 



lim w a {E) = —= — 



exp 



(E - E a - e a y 
2al 



where the quantities e a and a a are defined by 
e a = (a\H\a) - (a\H \a), 



al = (a\H 2 \a) 



(a\H\a) 



(7) 



(8) 
(9) 



Note that e a has a classical counterpart while a a is purely 
quantum mechanical. It appears because the commuta- 
tor [H, Ho] is nonzero, and the distribution w a (E) there- 
fore has nonzero width. 



The rigorous proof of equation Q is given in ^(J and 
is based on two assumptions: The energy of each group 
as defined in equation © has to be bounded, i. e. 



(x\nM < c 



(10) 



for all normalised states \x) and some arbitrary constant 
C, and 



a 2 a > N G C 



(11) 



for some constant C > 0. 

If these conditions are met, equation © can be com- 
puted for Ng > 1: 



(a\p\a) = 



exp 



(-0Va + W^ a ) erfc (Eo-y a + Pal 



2Z 



V2(Ja 



(12) 

where y a = E a + e a and crfc(a;) is the conjugate Gaussian 
error function. The second term appears because the en- 
ergy is bounded from below and the integration extends 
from Eo to oo. 

If the argument of the error function divided by V 'Ng 
is finite (different from zero), the conjugate error func- 
tion can be approximated by its asymptotic expansion 
for Ng ^> 1 erfc(x) ~ exp (—a; 2 ) / \/tt x for x — ► oo 
and erfc(x) «2 + exp (—a; 2 ) /V^x for x — > —oo. 

The off diagonal elements (a.|p|6) vanish for 
\E a — Eb\ > a a + a\j because the overlap of the 
two distributions of conditional probabilities becomes 
negligible. For \E a — Eb\ < a a + a^, the transformation 
involves an integral over frequencies and thus these 
terms are significantly smaller than the entries on the 
diagonal. 

We now test under what conditions the density matrix 
p may be approximated by a product of canonical den- 
sity matrices with temperature (3i oc for each subgroup 
p = 1,2,..., Ng- If we assume periodic boundary con- 
ditions, all reduced density matrices are equal and their 
product is of the form (a|/5|a) oc exp(— (3i oc E a ). Since the 
trace of a matrix is invariant under basis transformations, 
it is sufficient to verify the correct energy dependence of 
the density matrix in the product basis. Thus the loga- 
rithm of the rhs of equation (|12|) must be a linear function 
of E a plus a possible constant. It follows from the asymp- 
totic expansion of the conjugate error function, that this 
can only be the case for 



1 E a + e a -E - Pa 2 a 



V2a a 



>0, 



(13) 



meaning the lhs should have a finite positive value, and 



-£a + ^a~ C lK+C 2 , (14) 



where c\ and c 2 are constants. 

If |ci| <C 1 temperature becomes intensive: 

Ploc = d. 



(15) 
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Furthermore, to ensure, that the density matrix of each 
subgroup [i is approximately canonical, it must hold that 



l^wd^ + ca, (16) 



where e M = (a\V(q^ n , q^ n +i)\a), cr 2 = (a\Hl\a) - 

(a|W» 2 and u\ = T,£l_MH v -iHv + U v U v -x\a) - 
2{a\n v -i\a){a\H v \a). 

For a model obeying equations (|1(J|) and JTTJ, the two 
conditions l|13|) and (|16|) . which constitute the general 
result of this letter, must both be satisfied. These funda- 
mental criteria will now be applied to a concrete example. 

Thermal properties of insulating solids are quite suc- 
cessfully described by harmonic lattice models. We there- 
fore consider as an example a harmonic chain of N G ■ n 
sites with periodic boundary conditions. In this case, the 
potentials read 

u {<!%) = -J w o <$ and V(q it q i+ i) = -mw„ ft q i+1 , 

(17) 

where ft is the displacement of the particle at site i from 
its equilibrium position i ■ a with a being the distance 
between neighbouring particles at equilibrium. We de- 
vide the chain into N G groups of n particles each and 
thus get a partition of the type considered above. The 
Hamiltonian of one group is diagonalised by a Fourier 
transformation and reads in terms of creation and anihi- 
lation operators at and afc fl^ ]. 



Tin 



> ftm+1 ) 



y^t^fc ( a\a k + 
k 



(18) 



with k = Trl/(ao(n + 1)) (I = 1,2, ...,n) and 
cul = Alo 2 sin 2 (fca/2). We chose H = 1. 

We first check condition (|11|) . Expressing the group 
interaction V(g Mn , ftj n +i) also in terms of a k and a k , one 
sees that (7^ = for all fi and therefore a 2 = X)u=i a fi> 
where has a minimum value since [afe,aL = 6k P , 
(a 2 > Au> a n/2>Ti). a 2 a thus fulfills condition (JTTJ. 

To satisfy condition (|10|) we only consider states where 
the energy of every group, including the interactions with 
its neighbours, is bounded. Thus, our considerations do 
not apply to product states |o), for which all the energy 
was located in only one group or only a small number 
of groups. For N G ^> 1, the number of such states is 
vanishingly small compared to the number of all product 
states. 

To analyse conditions H13|) and 11611 , we make use of the 
continuum or Debye approximation [13J , requiring n>l 
and ao <C /, where I = nao, the length of the chain, is 
finite. In this case we have tok ~ v k with the constant 
velocity of sound v = ui a and a 2 = 4n~ 2 E^E^x- 

It is sufficient to fulfill conditions l|13fl an d i|16|) f° r 
some adequate energy range E m i„ < < E max . We 



therefore need to estimate E m i n and E max . The total 
chain is in a thermal state. Let E be the expectation 
value of its energy without the ground state energy, 



E = N G nk B <d 



e/T 



1 



dx. 



(19) 



Here, is the Debye temperature, which is a character- 
istic constant of the material chosen and can be found 
tabulated 01 . We then choose 



1 E 



< E u <a 
~ Ni 



E 

G 



— + ^ (20) 
N G 



with a 3> 1. The product of the strongly growing den- 
sity of states and the exponentially decaying occupation 
probability forms a sharp peak around E. The above en- 
ergy range has been chosen to be centered at this peak 
and to be significantly larger than its width. 

We now analyse conditions l|13|l and Ijl6|l for the har- 
monic chain in Debye approximation. In this case, con- 
dition (|13f) translates into 



E 



N G 



1 £ C 



E a -A[3n- 2 Y, 



i E^E^+i 



>0, 



y/2J4n 



where e M = 



E"G TP TP 

(21) 

dfj, = 0, n + 1 ~ n has been used. For 
E > Eq, condition i|16|) turns out to be stronger and l|21|) 
becomes irrelevant. For E < Eq, (|21|l is hardest to satisfy 
if all E^ are equal to the lower boundary of the energy 
range (|20J) . For n one thus gets 



n > 4 



6 a 
T i 



(22) 



where e = E / '(nN G ksQ) and we have used 
E = nN G k B Q/4. 

To investigate condition (|16fl we take the derivative 
with respect to on both sides, 



E^-i + E^+i — 2 



Nr, 



n 2 N G 



Cl 



(23) 



where we have seperated the energy dependent and the 
constant part on the lhs. 12-il) is satisfied if the energy 
dependent part, the first term, is much smaller than one. 
Taking £^_i and E^ +1 equal to the upper bound in equa- 
tion (|2nj), this yields 



2a 6 _ 



(24) 



where the "accuracy" parameter S <C 1 quantifies the 
value of the energy dependent part in J22Jl. 

Since, by virtue of equations Il_?2li and J22J, the con- 
stant part in the lhs of J53J) satisfies (2f3/n 2 )(E /N G ) < 
y/S/(y/2a) — S/a 2 <C 1, temperature is intensive. 
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FIG. 1: Log-log-plot of n m in from eq. i'2'It (dashed line) and 
Tlmin from eq. 1241 (solid line) for a = 10 and S — 0.01 as a 
function of T/Q for a harmonic chain. 8 and a are denned in 
equations 1211 and OH . respectively. 



Inserting equation l|19|) into equation (|22[l and 1|24|) one 

can now calculate the minimal n for given S, a, G and T. 
Figure shows n m - m for a = 10 and 6 = 0.01 given by 
criterion (1221) and (|24() . respectively, as a function of T/6. 



partitioned the chain into identical groups of n adjoining 
particles each. Taking the number of such groups to be 
very large and assuming the total system to be in a ther- 
mal state with temperature T we have found conditions 
(equations (|T3|> and ljll)|> ). which ensure that each group 
is approximately in a thermal state with the same tem- 
perature T, implying intensivity. In the quantum regime, 
these conditions depend on the temperature T, contrary 
to the classical case. 

We have then applied our general method to a har- 
monic chain. Equations (|25(l and (|26J) provide a first 
estimate of the minimal length scale on which intensive 
temperatures exist. This length scale should also con- 
strain the way one can meaningfully define temperature 
profiles in non-equilibrium scenarios flfjj . 

We thank J. Gemmer, M. Michel, H. Schmidt, M. Stoll- 
steimer, F. Tonner, M. Henrich and C. Kostoglou for 
fruitful discussions. 



For high (low) temperatures n r , 
mated by 



can thus be esti- 



2a/ 5 
(3a/27r 2 ) 



(9/T) 3 



for T > 9 
for T < 9 



(25) 



The characteristics of figure ^ are determined by the en- 
ergy dependence of a a - The low temperature behavior, 
in particular, is related to the fact that o~ a has a nonzero 
minimal value. This fact does not only appear in the har- 
monic chain but is a general feature of interacting quan- 
tum particles. The commutator [H, Hq] is nonzero and 
the ground state of the total system is energetically lower 
than the lowest product state, therefore a a is nonzero, 
even at zero temperature |a ll5j ■ 

Finally, the minimal length scale on which temperature 
exists according to our approach is given by 



In 



O , 



(26) 



where ao is the lattice constant. Temperature measure- 
ments with a higher resolution should no longer be inter- 
preted in a standard way. 

Let us conclude with some numerical estimates: Chos- 
ing the "accuracy parameters" to be a = 10 and 
S = 0.01, we get for hot iron (T > 6 w 470 K, 
ao ~ 2.5 A) l m in ~ 50 /im, while for carbon (9 ~ 2230 K, 
ao ~ 1.5 A) at room temperature (270 K) l m i n ~ 10 /im. 
The coarse-graining will experimentally be most relevant 
at very low temperatures, where l m i n may even become 
macroscopic. A pertinent example is silicon (O ~ 645 K, 
ao ~ 2.4 A), which has l m in ~ 10cm at T IK (again 
with a = 10 and 5 = 0.01). 

In summary we have considered a linear chain of par- 
ticles interacting with their nearest neighbours. We have 
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